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Abstract. We study aging dynamics in two non-disordered spin models with multi- 
spin interactions, following a sudden quench to low temperature. The models are 
relevant to the physics of supercooled liquids. Their low temperature dynamics 
resemble those of kinetically constrained models, and obey dynamical scaling, 
controlled by zero-temperature critical points. Dynamics in both models are thermally 
activated, resulting in multi-stage relaxation towards equilibrium. We study several 
two-time correlation and response functions. We find that equilibrium fluctuation- 
dissipation relations are generically not satisfied during the aging regime, but 
deviations from them are well described by fluctuation-dissipation ratios, as found 
numerically in supercooled liquids. These ratios are purely dynamic objects, containing 
information about the nature of relaxation in the models. They are non-universal, and 
can even be negative as a result of activated dynamics. Thus, effective temperatures 
are not well-defined in these models. 
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1. Aging dynamics and plaquette models 

When a liquid is cooled through its glass transition, its state changes: it becomes an 
amorphous solid yyj. In this non-equilibrium state, physical properties are not stationary, 
and the system displays aging behaviour IHj . A similar situation is encountered in 
many different materials, ranging from disordered magnets to dense granular media. 
A full understanding of the non-equilibrium glassy state remains a central theoretical 
challenge. In this work, we study the aging dynamics of two spin models with multi-spin 
interactions These finite-dimensional, non-disordered, dynamically heterogeneous 
systems have been of recent interest as models for the glass transition in supercooled 
liquids UUn\- 

Theoretical studies of mean-field models have given important insights into 
the aging dynamics of both structural and spin glasses [Hj. In these models, 
thermal equilibrium is never reached, and aging proceeds by downhill motion in an 
increasingly flat free energy landscape [H]. Time translational invariance is broken, 
and two-time correlation and response functions depend on both their arguments. 
The fluctuation-dissipation theorem (FDT), which relates equilibrium correlation and 
response functions, does not apply in the aging regime. Instead, correlation and response 
are related by a non-trivial fluctuation-dissipation ratio (FDR). This led to the idea 
that aging systems might be characterised by an effective temperature defined in 
terms of the FDR. Physically, relaxation in glassy systems occurs in well-separated time 
sectors; it is easy to imagine associating each sector with an effective temperature [TTj. A 
thermodynamic interpretation of effective temperatures has also been suggested |HIIT^. 
related to the concept of replica symmetry breaking. Taken together, these results lend 
considerable appeal to the the mean- field description of aging (see ^HI for a review). 

However, there are many systems of physical interest in which the dynamics are 
not of mean-field type, displaying both activated processes and spatial heterogeneity. 
Two examples are domain growth in disordered magnets and liquids quenched to 
below the glass transition ^Hl QB] • While some experiments and simulations [TS^ seem to 
detect a mean-field aging regime, theoretical studies have found ill-defined FDRs [T7\ . 
non-monotonic response functions dl dl 1201 HI]; observable dependence j22l 
non-trivial FDRs without thermodynamic transitions [Sj, and a subtle interplay between 
growing dynamical correlation length scales and FDT violations [211123; experiments 
have also detected anomalously large FDT violations associated with intermittent 
dynamics j2Sl- Moreover, at large times, aging often proceeds via thermal activation, 
and it was recently shown [22] that this can lead to negative response and well-defined, 
but negative, FDRs. 

In this work, our aim is to investigate further the effects of thermal activation and 
dynamic heterogeneity on aging dynamics. To this end, we study two different two- 
dimensional plaquette models with multi-spin interactions jH El El 1211 12H1 121] • On 
one hand, they can be viewed as finite dimensional, non- disordered versions of p-spin 
models, and can be viewed as an attempt to transfer mean-field concepts to the finite 



Fluctuation- dissipation relations in plaquette spin systems 



3 



dimensional world. On the other hand, they possess dual representations in terms of 
independent plaquette excitations with constrained dynamics, and are therefore directly 
related to the physics of kinetically constrained models (KCMs) jHO] , with the advantage 
that both spin and excitation degrees of freedom can be studied separately PT] . 

In these models, low temperature relaxation towards equilibrium proceeds via 
several distinct stages, each stage being associated with a particular energy barrier. 
While these separate time scales are superficially reminiscent of the time sectors found 
in mean-field spin glasses, their physical origin is quite different. We will show that 
the dynamics of plaquette models share some strong similarities with mean-field aging 
dynamics, but also important differences. In addition, the dynamics of these models 
becomes critical at low temperatures, where dynamical length scales diverge. Our work 
therefore also pertains to the study of FDRs in non-equilibrium critical dynamics [T3| I32j . 

The first model that we study is the triangular plaquette model (TPM) jHEHEHj, 
defined for Ising spins, = ±1, on the vertices of a triangular lattice. The Hamiltonian 
is 

-f^TPM = ^ S1S2S3, (1) 

1,2,3 e V 

where the sum is over the downward pointing triangular plaquettes of the lattice. We 
work with periodic boundary conditions in a system of linear size L of the form L = 2^. 
The dynamics involve single spin flips, with Glauber rates. 

It is useful to define the binary dual plaquette variables by 

ni = ^{l- S1S2S3), (2) 

where the three spins lie on the vertices of the ith downward-pointing triangular 
plaquette. The thermodynamics of the system are those of the non-interacting dual 
spins Hi. Combining (P) and ^ we have, at temperature T = the equilibrium 
defect density, (ni)eq = (1 + e^)~\ and an obvious connection to kinetically constrained 
models |11I2I]. 

The TPM has a hierarchical free energy landscape: as the temperature is lowered, 
the system falls into ever deeper potential wells, and the equilibrium relaxation time 
diverges as r ~ e^^^^"^^ [7j. This behaviour resembles that of the one-dimensional East 
model jHSl; or its generalization in higher dimensions The hierarchy of energy 

barriers in the TPM is clear from the behaviour on a quench: in Fig. ^ we plot the 
defect (or energy) density 

c(t) = (n.(t)), (3) 

as a function of time, after a quench from infinite temperature (the average is over initial 
conditions and over thermal histories). The energy barriers relevant for the quench take 
integer values, so the decay of the energy becomes a function of the scaling variable 
Tlogt, as long as t remains smaller than the equilibration time j36j . 
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Figure 1. Relaxation of the energy density in the TPM (left) and SPM (right), 
after a quench from infinite temperature to T = In the TPM, the relaxation 

proceeds in a hierarchical way, with stages corresponding to integer energy barriers 
[Hll^. The system has not yet equilibrated: this will occur when c{t) = (1 + e'^)"^. 
In the SPM, the decay has just two stages; the second plateau in the decay of c{t) 
represents equilibration of the system. The first stage of relaxation has no energy 
barriers, and corresponds to zero temperature dynamics (a quench to T — is shown 
in the inset, and compared with a simple exponential as a dashed line); the second 
step is characterised by activated reaction-diffusion behaviour 



Our second system is the square plaquette model (SPM) [2H1I11, in which we define 
Ising spins s, = ±1 on a square lattice, with 

H = ^ S1S2S3S4, (4) 

1,2,3,4 G □ 

where the sum is over the elementary plaquettes of the square lattice. For the SPM, we 
define the binary dual variables by 

Hi = ^(1 - S1S2S3S4), (5) 

where the four spins lie on the vertices of the ith square plaquette. We use the same 
symbol, rij, for the dual spins in both plaquette models. The thermodynamics are again 
those of non-interacting dual spins, but constraints from the boundary conditions cannot 
be avoided in this model [Zj. These are minimised by our use of periodic boundary 
conditions, but it is necessary to use linear system sizes L significantly greater than 
the inverse density of defects, L{ni) ^ 1. This is a rather strong constraint at low 
temperature. In the thermodynamic limit, the equilibrium defect density is the same as 
that of the TPM, {ni)^^ = (1 + e^)-\ 

The low temperature behaviour in the SPM is that of dilute point defects. In 
equilibrium, they diffuse at a rate that is proportional to e~^^; similar activated diffusion 
occurs in the Fredrickson- Andersen (FA) model [SZI. This process is mediated by pairs 
of defects that diffuse rapidly along one-dimensional paths that are aligned with the 
lattice axes. Where the TPM has a hierarchical relaxation towards equilibrium, the 
SPM has just two stages. The behaviour on a quench is shown in Fig. ^ The second 
stage is a reaction-diffusion process with activated diffusion, and the decay is a function 
of the rescaled time te~^^ [HI IE]- The early time regime of the SPM has no energy 
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barriers. It corresponds to relaxation into a jammed state at zero temperature. This 
relaxation is frustrated by entropic effects leading to a non-exponential relaxation of the 
energy density, as shown in the inset of Fig. [T] 

Earlier work on aging dynamics in the TPM and SPM |^ suggested that 
the combination of activated dynamics and hierarchical relaxation leads to novel and 
intriguing behaviour for response functions and FDRs. Here, we study these FDRs by 
computer simulations. We investigate spin and defect observables, including both local 
and spatially dependent correlations; we make use of a recently developed method for 
direct measurement of response functions in Monte Carlo simulations I39j . 

We find a rich structure in the FDRs, with qualitatively different behaviour in 
different stages of the relaxation. We observe FDRs close to unity for regimes in which 
energy conserving processes allow the system to explore configuration space efficiently; 
correlations may also relax by irreversible processes, whose rates are independent of 
the perturbing field; in this case the FDR is close to zero. The relative rates for these 
processes refiect local properties of the free energy landscapes in these systems jH]. We 
also observe negative FDRs, coming from activated aging behaviour, as in fI7\ . 

The rest of the paper is organized as follows: in Section II we study FDRs in 
the triangular plaquette model; in Section III we do the same for the square plaquette 
model; Section IV gives a brief summary of our results and their implications. Details 
of the no-field method for measuring response functions are provided in [Appendix A 



2. Fluctuation-dissipation relations in the triangular plaquette model 



2.1. Spin observables 

We begin by considering correlation and response functions for spins in the TPM. The 
two-time spin autocorrelation function is 

Cs{t,t^) = {si{t)si{t^)), (6) 

and its conjugate response is 

d 



(7) 



where hi is the strength of a magnetic field that acts on site i between times t^ and t. 
Averages are performed over realisations of the initial condition and the thermal history, 
as above. For these observables, the FDR, Xs(t,t„), is defined by 



dt^ 



-Xs{t, tv 



dCsit,t^ 



dt^ 



(8) 



In equilibrium, the fiuctuation-dissipation theorem states that X{t,t^) = 1 for all t and 
tw- However, away from equilibrium, there are no such restrictions, as we shall see. 

Following Refs. jHHl 1^11; two-time linear response functions can be computed 
without using a perturbing field (see [Appendix A for details). This allows direct access 
to the FDR. We emphasize that the derivatives in Eq. (|H)) are with respect to t^. For 
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Figure 2. (Top) Two-time spin autocorrelations (left), response functions (centre) 
and parametric FD plots (right) in the TPM at /? = 10. The correlation and response 
are plotted against t^, for various t. The FD plot (top right) uses as a parametric 
variable; the three traces at /3 = 10 are t = (70, 1.6 x 10^, 10*). The equilibrium FD 
relation, xititw) = 1 — C{t,tvj), is shown as a dashed hue. We also show parametric 
plots for /3 11 and t = (70, 1.7 x lO"^). Since these values of t are within plateaux 
of c{t), the FDR depends very weakly on temperature, although the correlation and 
response have quite different time dependence (not shown). (Bottom) Sketches of these 
plots, illustrating the various definitions given in the text. Note that the parameters 
Xnm label the gradient of the parametric plot. 



this reason, we present data for fixed time t, as a function of t^, unlike previous work 
on the same models where the opposite convention was used [3 E] . We shall see that 
this difference qualitatively affects the results and their interpretation. 

The dependence of the FDR on t„ can be replaced by a parametric dependence on 
the value of the correlation: 

^(9.^) = ^(^.UIc(Mw)=.- (9) 
That is, working at fixed time t, the waiting time t^, is parametrised by the value of 
C(t, tw)- A parametric plot of xil, t) versus q (an "FD plot") has local gradient —X{q, t), 
corresponding to the FDR defined in Eq. (jH)). FD plots were suggested by the study of 
mean-field models where Eq. simplifies to a single argument function X{q, t) = X{q). 
This asymptotic property has not been confirmed beyond mean-field |T^ . 

The behaviour of the two-time spin autocorrelation function of Eq. ^ is shown in 
Fig. 121 If t is held constant, and t^ decreased, the correlation decreases, in stages that 
mirror the decay of the energy after the quench (recall Fig. P). The generic structure of 
the correlation is sketched in Fig. |21 

The plateaux in the energy, c(t) , come from the system being trapped in metastable 
states. The dynamics are effective in exploring the metastable state: in the nth such 
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state, events corresponding to energy barriers smaller than n are very fast, while those 
corresponding to energy barriers greater than or equal to n are frozen, until such times 
that the system escapes from the metastable state. Thus, the autocorrelation function 
quickly relaxes to a value g„„, which measures the self-overlap of this metastable state. 
As tw is decreased further, we arrive at a new plateau in the the two-time function. This 
measures the mutual overlap between the plateau state n and the previous plateau state 
n — 1. We denote this quantity by g„,n-i; earlier plateaux in the two-time function give 
the mutual overlaps g„m with m < n, see Fig. |21 

The two-time local integrated response, Xs(t,tw), is also plotted in Fig. |21 The 
stages of the dynamics are still visible, although the response decreases with increasing 
t^, unlike the autocorrelation function. Note also that curves for different t cross one 
another at short t^. This means that response functions plotted at fixed t^ for increasing 
t are non-monotonic, which illustrates why the use of t as a parametric variable in FD 
plots produces quite unusual results [T71I2I|. We sketch the behaviour of the response 
function in Fig. El and define the plateau values of the response function by Xnm- By 
analogy with the mutual overlap, the mutual response Xnm is the value of the response 
function measured in the stage n, given that the field was continuously applied from 
stage m to stage n. 

If the response is now plotted parametrically against the correlation, then each 
stage of the dynamics appears as a separate segment on the FD plot. It is clear from 
Fig. 121 that X{q,t) 7^ 1 in the TPM, reflecting the fact that the system is well out of 
equilibrium. More importantly, it is apparent that the segments of the parametric plot 
in that figure are close to straight lines, so that FD plots are reducible to a discrete 
set of numbers. We denote the slopes observed in the FD plots by Xnm, as sketched in 
Fig. El We define 

J- Qnn 

and 

Xnm Xn,m+1 (T -\\ 

for n > m. These equations are a discretized version of Eq. i^. Physically, each stage 
n of the hierarchy corresponds to a given timescale, t„ ~ e"/^, and a well-defined length 
scale, in, giving rise to well-defined two-time (or rather 'two-stage') correlation and 
response functions and FDRs. The whole dynamical behaviour is therefore encoded 
in a discrete set of real numbers, {qnm,Xnm,Xnm}, which obey certain relations, as we 
shall see in section 12. HI The FD plots found for spin observables have similar shapes to 
the ones reported in Ref. j2ZI for the East model, but differ significantly from earlier 
work ini- 

The first segment of the response always has a gradient very close to unity, x„„ ~ 1, 
consistent with quasi-ergodic (or quasi-equilibrium) dynamics within the metastable 
plateau state. We also note that the last plateau, which corresponds to correlation and 
response between stages n and 0, has x„o — 0. This is to be expected, since the first stage 
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of the quench only involves unactivated relaxation events, and these are accepted with a 
probability close to unity. Thus, their rate depends very weakly on any perturbing field, 
and the response associated with this completely irreversible dynamics is very weak 
(more precisely, it is 0{e~^), which does indeed vanish at small temperature). This is 
very similar to the zero-FDR found in coarsening ferromagnets I41j. 

We have discussed the data in Fig. |21at a single temperature. However, the scaling 
behaviour of the TPM with temperature is well- understood [71 We find that the 
plateau values of the energy, the self-overlaps qnm and the ratios Xnm depend only very 
weakly on temperature. As the temperature is reduced towards zero, the stages of the 
dynamics become more clearly defined and timescales more separated. If we plot data 
as a function of the rescaled times, (z/, z/^) = (Tlogt, Tlogt„), then the plateaux fall on 
top of one another, and the parametric plots quickly become independent of temperature 
(for a given value of z/). This assertion is confirmed numerically in Fig. |21 where we 
show that parametric plots for different temperatures, /3 = 10 and f3 = 11, but similar 
u, perfectly superimpose. 

Mean-field studies have suggested the possibihty to define effective temperatures 
from FDRs through 

= (12) 

Our finding of piecewise linear FD plots in Fig. |21 apparently offers a nice illustration 
of the intuitive idea that well-separated relaxation timescales could lead to quasi- 
thermalization at a given effective temperature within each time sector, as is the case in 
mean-field spin glasses jHI. A further property derived from mean-field models is that 
the value of Tgs is shared by all physical observables, physically implying that different 
degrees of freedom have also thermalized between themselves. We now investigate this 
issue. 

In the TPM, the spatial dependence of two-spin correlation functions is trivial, 
since both two-spin correlations and responses are purely local. The symmetries of the 
model imply |3] 

{si{t)s,{t^)) = 6,,C,{t,t^), (13) 

and 

^{sj) = 5ijXsit,t^), (14) 

which trivially implies that all two-spin correlation functions are equivalent. 

In order to study spatial correlations, it would in fact be necessary to consider 
four-point functions the relevant correlations were discussed in j7], but we restrict 
ourselves to two-point functions in this paper. Interestingly, ()14|) also implies that the 
dynamics of the total magnetization are trivially related to those of local observables. 
This remark is potentially relevant for supercooled liquids, where static two-point 
correlations are also believed to be decoupled from the microscopic physics of cooperative 
motion. 
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Figure 3. Parametric FD plot for defects in the TPM, at /? = 10 and t S 
{70, 1.7 X 10^, 10^}, with t decreasing from top to bottom. The structure is the same 
as for the spin observables in Fig. |21 but the FDRs are quantitatively different. The 
dashed hne is the FDT: Xd(i,^w) = 1 — Cd(i,^w)- 



2.2. Defect observables 

To understand the significance of tlie FDRs (and tfierefore effective temperatures) 
revealed by spin observables, it is natural to ask whether different physical quantities 
exhibit a similar behaviour. We therefore consider the dynamics of the dual plaquette 
variables defined in Eq. ((21) . An excited plaquette, or "defect", is a triangular plaquette 
for which Ui = 1. Those defects become increasingly dilute at low temperatures. We 
consider the defect autocorrelation function, 

^ / , , X _ (n^(t)n^(tw)) - {njjt)) {ni{t^)) 

(n.(t))[l-(n,(t))] ' ^'"^ 



and its conjugate response, 
Xd{t,t^) 



dm 



■{n^{t)) 



(16) 



where the perturbation SH = —QiUi acts between t^ and t, and we evaluate the right 
hand side at gi = 0. We use connected correlation functions and we have normalised both 
correlation and response, so that Cd(t,t) = 1, Cd(t oo,t^) = 0, and the fluctuation- 
dissipation theorem for equilibrated systems reads Xd(t,t^) = 1 — Cd(t,t^). 

Overall, the results for defects are very similar to those for spins: dynamic functions 
can be reduced to discretized overlaps and responses. Moreover we find FD plots made 
up of linear segments, as shown in Fig. IHl However, we find that the FDRs associated 
with the defects differ quantitatively from those of the spins, so we conclude that these 
numbers do not lead to useful effective temperatures. 

These results are similar to those found for the East model |2Z1. However, the 
low temperature hmits of the g„m and Xnm can be estimated in the East model by 
rather simple arguments whereas those in the TPM are more complex, and contain 
information about the mechanisms of relaxation in that system, as we now discuss. In 
the East model, relaxation from stage n — 1 to stage n effectively consist of selecting a 
random set of excitations and removing them. This can be represented schematically 
by: 
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Figure 4. Processes in ttie TPM requiring an energy barrier of 1 unit. Spins are 
defined on the vertices of the lattice. The x signs mark defects. The black dots 
show spins that must flip for a transition between the left and right states. (Top) 
Irreversible coagulation: a downward direction in the energy landscape. (Bottom) 
Reversible ('diffusive') event: the energy landscape is locally flat. 



1. .1.1.1.1.1. .1.. 
1..1...1...1..1.. 

1 1 1.. 

where we show configurations of the East model's defect variables, at three times, with 
the latest time at the bottom. Neighbouring defects in the first state are separated by 
gaps of length at least 1, so it representative of the first plateau. The second and third 
states are representative of the 2nd and 3rd plateaux respectively (gaps at least 2 and 
4). We can see that 

{ni{t)ni{t„)) (mit)), (17) 

so that 

„(East) ^ max ^ 1 - (^»(^w)) 

'dnm -'dnm - ^ _ ^^^^^^^ , I ) 

where t is a time within the nth plateau of c{t), and tw a time within the mth plateau. 

In the TPM, the situation is different, since there are both reversible and irreversible 
moves taking place during the relaxation. Irreversible coagulation events involve 
removing two defects and adding a new one on a different site. Schematically, relaxation 
proceeds as: 

1..1..1..1.1..1. 
1...1....1...1.. 

where pairs of defects combine to leave a single defect somewhere between the previous 
two. A more realistic representation of one such process is given in the top panel of Fig.HJ 
Additionally, reversible moves such as the one shown in the bottom panel of Fig. H] can 
occur. In that particular example, two defects separated by distance 2 'rotate' around 
a given plaquette. This process occurs on the same timescale as the coagulation event 
shown in the same figure. Both types of moves reduce the self and mutual overlaps of 
the plateau states. Numerically we confirm that q^™ < qnm^\ but we are not able to 
estimate analytically. 
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The second difference between the TPM and East models is that dynamics at low 
temperatures in the East model are those of irreversible events which allow no response, 
so that defect FDRs vanish as T — llHj. In the TPM, the Xnm are finite numbers 
(except Xno — 0); they become independent of temperature as T — 0. The rate for 
'coagulation' events is independent of the local field at low temperatures, but 'diffusive' 
events do allow a response. For instance, if we take the top left state of Fig. HI as our 
initial condition, then we end in the top right one regardless of the local field: the 
response is zero. On the other hand, if we start in the bottom left state then the final 
state may be either the bottom left or bottom right one, with probabilities weighted 
by the local fields: the response is finite. In the energy landscape picture F^], 'diffusive' 
processes are along directions in which the energy landscape is fiat, and 'coagulative' 
processes correspond to downward paths in the energy landscape. We conclude that 
Xn,n-i measures the relative rates for 'diffusive' and 'coagulative' events during the 
relaxation: a large value of Xn,n-i is a sign that 'diffusive' processes are dominating. We 
will confirm this statement explicitly for the SPM in section El since the relative rates 
for the two types of process can be controlled in that model. 

Finally, we note that the field Qi is a local perturbation to the temperature; in the 
East model, this leads to negative FDRs, if the time t is between plateaux of c(t) p7| liH]. 
The same behaviour occurs in the TPM, but we consider only the plateau FDRs, Xnm, 
in this article. These quantities are positive (except for x„o — 0). 

2.3. Independent steps approximation 

In the sketches of Fig. |2l we reduced the behaviour of local dynamic observables to a 
discrete set of numbers. We now show that if successive stages of the dynamics are 
statistically uncorrelated, then several relations between these numbers can be derived. 
These relations are satisfied by our numerical data. 

We use only the fact that probabilities for combinations of independent events 
factorise. We define Ps{t,t') to be the probabihty that a defect survives on a given site 
from time t' to time t > t'. If relaxation proceeds by independent stages then we expect 

Psitnjtm) = Y\. Ps(tn',tn'-l), (19) 
n>n'>m 

where the times {ti} separate the stages. We consider a situation where t„ is a time 
within the nth plateau. 

A naive approximation for the survival probability would be to assume that 
Ps(t, tw) — {nunit^) / {uit^) , which is the probability that there is a defect on site 
i at time t, given that there was one on site i at time t„. In the East model, 
this approximation is appropriate: we know that qnm = 'y^rrT) it follows that 

{uit^nuj ~ {uitj nn'=m+iW„'^iV+i)/(^«v)- 0*^6^ ^a^d' the approximation 

fails if successive plateau states are completely uncorrelated, in which case {nitUu^) = 
{nit){nit„). In this latter case, the naive approximation fails because it includes in the 
survival probability events in which the defect on site i is destroyed, and then replaced 
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Table 1. Mutual overlaps and FDRs in the TPM at /3 = 10. We test the predictions 
of Eqs. H22|) and (|25|) . Consistent with those equations, we find q^i ~ (j'32'?2i and 
X21 — X31. We estimate the uncertainty for the spin overlaps to be less than 10% and 
those for the defect overlaps to be less than 20%; those for the FDRs are given. The 
(*) indicates that the defect value of 2:31 was evaluated at t — 10* (as in figure |2Jl; 
the third stage of relaxation was not complete at this time, but the FDR still seems 
consistent with 



921 



932 



931 



(932921 



X21 



X31 



Spins 
Defects 



0.63 0.43 0.26 (0.27) 
0.25 0.053 0.009 (0.013) 



0.25 ±0.02 0.21 ±0.04 
0.18 ±0.02 0.17±0.03(*) 



by a new defect. To take this into account, we estimate instead 

p / , , \ {njtUit^) — {nit){nitj) 

(-.*J(1-KJ) • ^ ^ 

With this approximation, the factorisation of probabihties for independent stages leads 
to 

Cd{t, tw) - Cd{t, t„_i)Cd(t„_i, tw), (21) 

which follows from p9|) and (f^UI) . if the dynamics between times t^-i and t are 
statistically uncorrelated with the dynamics between times t^ and Times t and t^ 
are chosen within the nth and mth stages, respectively. Correlators for spin observables 
factorise in a similar way. In terms of overlaps between plateau states, it follows that 

<inm = Y\. ^n',n'-l- (22) 

n>n'>m 

This relation is satisfied both for the East model and for completely uncorrelated stages. 
The dynamics of the TPM are rather slow, so that our data is limited, but we do obtain 
numerically that ^31 ~ 532^21 for both the spin and defect autocorrelations, see Table [U 
To establish the effect of the independent stage approximation on the FDR, we 
consider the impulse response. In the presence of an instantaneous field h acting at 
time t^, we expect the response to be given by the density of field-induced defects at 
stage n — 1, multiplied by the probability of these defects surviving to time t. That is, 
c%iW)?»_p/, , ^ d(n^(t„_i))fe 

~^ir~ - ' (23) 

which holds if relaxation between and t is independent of that between t„ and tn-i, 
as before. The derivative of x(t, t„) with respect to t„ is a response to an instantaneous 
field at ty^: taking care of normalisation, we find that 



dXd{t,t^ 



dt. 



(24) 



dt^ 

Again, the analysis for spin observables is similar. Using the definition of the FDR, and 
Eq. (PT|) . we arrive at X{t,t^) = X(tn-i,t^), and hence 

•^nm -^n— 1,771) (25) 
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for n — 1 > m. Thus, Xnm is a function of m only. The 'self FDR', Xmm, is close to 
unity. Physically, the FDR only depends on the stage of the dynamics during which 
the instantaneous field is applied. Our numerical results are again consistent with this 
analysis, see Table [T] 



2.4- Spatial structure in defect FDRs 

In contrast to the spin degrees of freedom, defects have non-trivial spatial correlations 
in the TPM (except at equilibrium, where defect correlations vanish at equal times). 
Defining the Fourier transform of the plaquette field nk{t) = N~^^'^ ?T,j(t)e*'° '^% where 
N is the number of plaquettes, we consider the correlator 

~ E|fc'|<fc(^fc'(^)^-fc'(^w)) - Sk'jrik')^ 

Gc((fc, r, r„j — T— 

Z^|fc'|<fcUJ 

= ^ E^'^d^^ - ^^-1) [(^^W^.(^w)) - (nMMt.))] , (26) 

where Ffc(|r|) = ^''^'^ ^^-^^ decays on a length scale of the order of from a value 
of unity at the origin. (Sums over wave vector k are over the first Brillouin zone, with 
additional restrictions as specified.) Compared to the apparently simpler correlator 
(nfc(t)n_fc(tw))5 Cd{k,t,t„) attaches more weight to short distances; measurements of 
this quantity are therefore less affected by noise in the correlations at large distances. 
Moreover, Cd{k,t,t^) interpolates smoothly between the local autocorrelation function 
when k is large, and the global correlation function of energy fluctuations at A; = 0. 
Normalising as in the previous sections, we define 

Crf(fc,t,0 = SiMlM (27) 
Cd{k,t,t) 

We begin with the case t = t„, which defines the static structure factor for the 
defects at time t^j- 

S{k,t^) = Cd{k,t^,t^). (28) 

We plot this function in Fig. El We see that S{k, t^) is an increasing function of /c, 
which means that defects are surrounded by regions of reduced defect density. This 
is consistent with a local coagulative process which induces some effective repulsion 
between defects. There is a length scale associated with the fc-dependence of S{k,t^): 
the time scale associated with coagulation of defects increases with their separation, so 
this length scale increases with t^- 

Moving to the time dependence of the correlation functions, we compare correlators 
at different wave vectors in the middle panel of Fig. |S| The plateau structure that is 
seen in the local correlations is not present at small k. Instead, the self overlaps of 
the plateau states approach unity, and the mutual overlaps approach zero. This is to 
be expected: within the plateau state, relaxation happens on small length scales, so 
fluctuations on large length scales do not relax, and self overlaps are large, for small k. 
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Figure 5. (Left) Structure factor of the defects, S'(fc,i„), at (3 — 10, and with t^, 
logarithmicaUy spaced between 10~^ and 10^ (from top to bottom). A growing length 
scale is apparent, whose origin is the increasing range of coagulative events at large t^. 
(Center) We show the correlation function of Eq. for k e {fcmax, fcmax/10}, where 
kmax is the largest wave vector in the Brillouin zone. We have lowered the temperature 
to /3 = 20, and restrict ourselves to the first and second stages of relaxation. (Right) 
FD plots at /3 = 20 for k G {fcmax, 2fcmax/3, fcmax/2} (from top to bottom) showing 
decreasing response with decreasing k. The time is i = 2 x 10^°; at this temperature 
then this time is within the second plateau of c{t). 



On the other hand, successive plateau states do not retain any memory of the average 
density in the preceding plateau, so mutual overlaps at small k are very small. 
The normalised response conjugate to Cd{k, t, t„) is 

It measures the response to a random field which is correlated over a length scale of 
k~^. We present FD plots for three wave vectors in Fig. at a single time t that is 
within the second plateau of c(t). We find that the amplitude of the response functions 
quickly decrease when k decreases and vanish as A; — > 0. This can be simply inferred 
from Fig. [T] since the energy density of each plateau states is temperature independent. 
Hence, the linear response of the energy to a temperature change vanishes in the small 
temperature limit. The situation is therefore similar to the one found in Ising or FA 
models [SHI /c-dependent FD plots smoothly interpolate from the one obtained 
from local autocorrelation functions to the one obtained for global quantities at A; = 0, 
the crossover taking place when k^(t^) ^ 1, where ^(t„) is a typical length scale 
characterizing the correlations between defects in the system. 

3. Fluctuation-dissipation relations in the square plaquette model 

3.1. Early time regime: zero temperature dynamics 

We now turn to the SPM, in which there are two relaxation stages. We begin our 
discussion with the initial stage, shown in the inset of Fig. ^ The corresponding FD 
plots for local spin and defect dynamic functions are shown in Fig. IHl 

For the defect variables, we find that the self overlap within the plateau state 
is gii ~ 0.82. This number comes from isolated spins that can fiip without energy 
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Figure 6. FD plots for the first stage of relaxation in the SPM at /3 = 20. We show 
the FDR for t = 5000, which is within the plateau of c{t), and for several earlier times, 
at which the system is approaching the plateau. (Left) Defects, (right) spins. We show 
the equilibrium FDT with a dashed line. 



penalty [Sj. The mutual overlap between the plateau and initial states is very close to 
zero, gio — 0.01. A well-defined FDR can be determined, Xio — 0.35, see Fig. El This 
makes the FD plot apparently very similar to the one found during the aging dynamics of 
mean-field spin glasses characterized by a one-step replica symmetry breaking transition 
at the static level [3 El El- 

For the spin degrees of freedom, we find qu ~ 0.98 and xn ~ 1. For smaller 
correlation, however, the parametric FD plot is curved and there is no well-defined 
value for xio- 

This behaviour is qualitatively similar to that of the FA model, characterized by 
an initial "zero temperature" relaxation, followed by a reaction-diffusion aging regime. 
However, there are some important differences. Firstly, the decay into the initial plateau 
is exponential in the FA model, and the overlap Qiq takes its maximal possible value 
gio ~ , _ — , where the defect density in the plateau state, and no = 1/2 

J- ^plateau ^ 

is the density in the initial state. This reflects the situation described by Eq. (fTH|). in 
which the decay proceeds by uncorrelated relaxation events. Similarly, the FA model 
has xiQ — > at low temperatures: this occurs because all Monte Carlo acceptance 
probabilities are either zero or unity at T = 0, and the dynamics are independent of the 
perturbing fleld; thus there is no response. On the other hand, the SPM has a flnite 
value for xio even at zero temperature, because there are energy conserving processes 
that have fleld-dependent acceptance probabilities close to 1/2, even when T = 0. 

We believe that the non-trivial zero temperature FDRs for both spins and defects 
come from the interplay between two types of process: reversible 'diffusive' spin flips in 
which the energy does not change, and irreversible 'coagulative' spin flips in which the 
energy is reduced. Both processes are occurring with similar rates. As discussed for the 
TPM, if the reversible process was dominating we would expect a large value for xio, 
while irreversible relaxation would lead to a small value of Xio- 

To confirm this analysis and probe in more detail the competition between diffusive 
and coagulative processes, we modify the Glauber dynamics, by introducing the 
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Figure 7. Correlation and response at T = 0, and their dependence on 72. 
(Left) Energy density c{t) (linear scale) after a quench to T = 0. Two stages are 
clearly visible at small 72; at large 72, the extra energy-conserving processes allow 
more efficient exploration of configuration space, and the energy relaxes to a smaller 
value. (Right) Parametric plot for spin observables at the latest values of t. At small 
72, the first stage of relaxation is purely 'coagulative' so it has x ~ 0; the second stage 
has both 'diffusive' and 'coagulative' processes, so x is finite. At large 72, energy- 
conserving (diffusive) processes dominate, and the FDR is large. 



additional parameters 7 in the rates of Eq. (jA.ljl (see [Appendix A ) . Our aim is to 
affect the dynamics of the SPM in such a way that diffusive and coagulative processes 
become well-separated processes, in order to analyze the consequences for the resulting 
FD plots. 

There are five types of move in the SPM: the energy change may be any of ±2, ±1, 0. 
We choose 7 to depend on the modulus of the energy change and label these rates 
according to the number of defects adjacent to the spin that flips, u. This number 
satisfies Wfinai = 4 — Winitiai, SO we define three 7ui„itiai factors for the three possible cases: 

(■u = 4) ^ (m = 0) multiplier 74; 
(m = 3) ^ (m = 1) multiplier 73; 
{u = 2) ^ {u = 2) multiplier 72. 

The multipliers are the same for both forward and reverse processes, since the dynamics 
respect detailed balance. At zero temperature, the rates for the three forward processes 
are all independent and are given by (74, 73, (72/2)). We find that the behaviour depends 
only weakly on 74 as long as 73 > 74. We therefore set 73 = 74 = 1, so that 72 is a 
dimensionless measure of the rate for the energy- conserving, diffusive, events. 

We show the energy decay and corresponding FD plots for spin observables in 
Fig. m as a function of 72. For 72 <^ 1, the zero temperature relaxation takes place in 
two stages, unlike the case of Fig.Hl in which 72 = 1. The first stage corresponds to pure 
coagulation: the system takes the path of steepest descent in the energy landscape. This 
process has an FDR of zero, as was the case for the first stage of relaxation in the TPM. 
In the second stage, the SPM explores its energy landscape by decreasing the energy 
where possible; otherwise, it makes energy-conserving moves. This stage has a finite 
FDR that reflects the ratio of energy-conserving, diffusive events (along flat directions 
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of the energy landscape) and coagulative events (along downward directions). We also 
note that this second stage is not a simple exponential decay, so the structure of Fig. IHl 
is not simply a result of two exponential relaxation mechanisms. 

In the case 72 ^ 1, the increased rate for diffusive moves means that the system 
prefers to explore the energy landscape in directions along which the landscape is flat. 
It therefore explores more of configuration space, and is more likely to find pathways to 
low energy. Thus, the plateau state contains fewer defects than that for 72 = 1. Since 
the diffusive moves explore the constant energy surface in an unbiased way, the system 
tends to equilibrate locally between each coagulative move, and the FDR approaches 
unity as 72 is increased. 

Overall, we see that the finite value of Xio observed in Fig.lHlcomes from competition 
between diffusive and coagulative processes. Its value, xio — 0.35 for 72 = 1 in fact 
depends continuously on the microscopic rates in the problem: we conclude that it is 
a non-universal number that measures the extent to which coagulation dominates over 
diffusion in the relaxational dynamics. In particular, the FDR has no connection with 
static quantities, in contrast to mean- field models [HI [12] ■ 

Moving to spatial structure, Eqs. ([T!?|IT^ also hold for the SPM, and spin correlation 
and response functions are independent of wave vector k. In the short time regime, the 
defect response is 0{k) at small k. This is because the zero temperature response comes 
only from energy-conserving processes, and the rates for these processes couple only to 
temperature gradients, and not to the absolute temperature. Our numerical simulations 
confirm that the response vanishes at small k. This was already demonstrated for the 
TPM in Fig. [S] so we do not show more data for the SPM. 

3.2. Aging regime: activated dynamics 

When the SPM is quenched to a small but finite temperature, the system eventually 
leaves the jammed state into which it relaxed initially, and enters an aging regime, in 
which the energy decays as a power law. We show local correlation-response data for 
both spins and defects in Fig. [HI As t and t^ get large, the local FDR approaches unity 
for all values of C, even though the system is still out of equilibrium. 

In the framework of the previous sections, the approach of the FDR to unity comes 
from the fact that defects are very sparse at large t^, and coagulation events become 
increasingly rare. The result is that diffusive processes dominate the response, and lead 
to X{q,t) ~ 1, for finite q and large t [recall that X{q,t) is the FDR, evaluated at 
C{t, tw) = 9; see ©]. However, it is now a familiar feature of react ion- diffusion systems 
that there is structure in X{q,t) for finite t and small q [23 [211 • define 

Xoo(t) = limX(g,t). (30) 

It was recently shown [2Z| that in the aging regime of the FA model, Xoo{t) is a negative 
number that depends on the dimension of the system, but not on the time t. Observables 
associated with domain walls in the Glauber-Ising chain have similar behaviour, with 
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Figure 8. Correlation and response in the aging regime for the SPM. (Left) Defects: 
we show t e {140, 2500, 1.0 x 10*, 4.3 x 10^, 1.8 x 10^ 7.4 x 10^}, which decrease from top 
to bottom. (Right) Spins: final times are t G {140, 1.6 x 10^, 1.1 x 10^, 7.4 x 10^}, with 
the same trend. For both cases, the first time is representative of the initial plateau in 
c{t): recall Fig.|Bl Except for the latest time, we concentrate on the region in which 
equilibrium FDT is not obeyed; the parametric plot goes through (C, x) = (IjO); for 
all values of t, but we do not plot data in this region. 



Xoo = jSHj. However, for both these models, X{q,t) only converges to its limiting 
value at very small q. This requirement also necessitates working at large t: the limit 
of long times and very small q makes direct investigation of this limit very difficult in 
simulations. For the FA model and the Glauber-Ising chain, Xoo is independent of the 
observation wave vector, k. From the point of view of simulations, it is fortunate that 
Xoo is typically much easier to measure at small wave vector. 

In Fig. El we show that the spatial structure of defect correlation and response in 
the SPM is consistent with this behaviour. That is, X{q, t) is positive and close to unity 
for large q, and negative for small q. The value of q at which the crossover takes place 
depends on the observation wave vector k; its value is close to unity at large k and close 
to zero at small k. The physical origin of a negative asymptotic FDR is the same as for 
the FA and East models j2Z]: when applying a small temperature change 6T > 0, (T is 
the field conjugate to the energy), the dynamics of the system gets accelerated and the 
energy decays faster towards equihbrium so that Sc(t) < 0, and indeed Sc(t)/ST < 0, 
unlike the equilibrium case. 

We now discuss how the limiting negative value of the FDR and the overall 
behaviour observed in fig.lHlcan be predicted by a similar calculation to that used for the 
FA model in Ref . |2Zj • The calculation for the SPM is slightly more involved, due to the 
presence of two growing length scales in the aging regime of that model. As discussed in 
in], the reaction-diffusion dynamics involve exchange of defect dimers between isolated 
single defects. This process takes place over a length scale i{t) ~ [l/uiit)]; at long times 
then this length scale is much larger than the correlation length for density fluctuations 
in the aging regime, which is ^(t) ~ [l/ni(t)]i/2. The processes by which dimers are 
exchanged have an energy barrier of two. Their rates are also inversely proportional to 
the distance i{t) over which the dimer must travel 01^1 ■ The result is that these rates 
depend linearly on the defect density in the system, averaged over a region of linear 
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Figure 9. Parametric FD plot for spatially dependent defect observables at /? = 5 
and t = 7.4 x 10^ (this time is within the aging regime). The wave vectors are 
fc/fcmax e {1-0, 0.50, 0.33, 0.24, 0.20, 0.17, 0.14, 0.11, 0.083, 0.050, 0.0}, with k decreasing 
from top to bottom. Dashed lines show the equilibrium FD relation, and the small k 
prediction of (jB.29|l . The deviations from the small k prediction arise from corrections 
to the aging limit, as discussed in the main text and in Fig. 1101 The response at fc = 
was calculated using an explicit field, rather than the no-field method described in 
[Appendix A| since that method is rather inefficient for responses at small k. 



extent C.{t). 

In [Appendix B[ we show how this situation can be described by an effective reaction- 
diffusion system, where correlations in the aging regime can be calculated in a field- 
theoretic formalism. The rates for reaction and diffusion depend the density, averaged 
over a region of fixed size I. We neglect the time dependence of this length scale, but 
our results depend only on whether k^ is large or small compared to unity. Hence, our 
results apply if ki(t^) and ki{t) are both small, or both large. The other parameter 
that enters the FDR is Zk{t^) = Dk'^/[2Xn{t^)], which measures the relative size of k 
and ^(t). 

The calculation described in [Appendix B is quite lengthy, but the important result 
is ()B.29j) . At small k, the FD plot is a straight line with gradient X = — (5/2). At 
large k, we expect quasi-equilibrium behaviour, X = 1 (reversible 'diffusive' processes 
dominate on small length scales). Further, we can identify an intermediate regime, in 
which ki(t) and k£{t^) are both large, but k^(t) and k^(t^) are both small. We find 
that this regime has X = —3. 

Comparing with the results of figure IHl the behaviour is in qualitative agreement 
with the field theoretic calculation. However, quantitative agreement is rather poor 
for the limit of small wave vector, and the intermediate regime is not apparent. To 
understand this, we observe that the field theoretic calculation is valid only deep in the 
aging regime, which requires 

(n,(t„)), imit)) < 1, imit)) > e-f^ (31) 

We also predict logarithmic corrections to these tree level results, since the critical 
dimension of the reaction-diffusion system is = 2 (to show this, one must be careful 
to take n(t)i(ty to infinity, in the limit of large time). 
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Figure 10. Parametric FD plots in the SPM at fc = 0. As the bounds of (|31|l are 

improved, the data approaches the prediction X = —(5/2), which is shown as a dashed 
line. (Left) Increasing (3 with a constant value of the rescaled time (c^t = 3.3) improves 
the bound {n^(t)) > e"'^. This improves the agreement with theory for t„ close to t. 
(Right) Increasing t and tw at fixed temperature improves the bound (ni(t„)) <C 1, 
which increases the range over which the system is close to the prediction for the aging 
limit. 



The conditions of (jM]) are rather restrictive, necessitating simulations at low 
temperatures and long times. While we were not able to access the regime of 
intermediate k, we show in figure^! that the FDR at A; = does appear to approach 
— (5/2) as we improve the bounds of (j!TT|) . Our data is therefore consistent with a value 
of = —(5/2), although convergence to this limit is rather slow. 

4. Conclusion 

Finally, we bring together the results of the previous sections and discuss their 
significance. We have shown that the FDR in these plaquette models has considerable 
structure. By using the 'no- field' method of Refs. [3111311, we have been able to perform 
systematic measurements of correlation-response relations allowing resolution of the 
questions about the SPM and TPM that were raised in Refs. [51 121j. 

By definition, the FDR measures the ability of a system to respond to an external 
perturbation (normalised by an equilibrium system with the same correlations). We have 
argued throughout that when the temperature is small, the system can only respond 
by 'diffusive' energy- conserving processes. A large FDR means that these processes are 
dominating (for example, this occurs in dilute reaction-diffusion systems); a small FDR 
means that irreversible 'coagulative' processes are dominating (for example, the East 
model). This is consistent with our interpretation of the FDR as a purely dynamic 
object in these systems, which have a single pure state at all finite temperatures. It 
contrasts with the situation following a quench into a regime in which ergodicity is 
broken, where the long time limit of the FDR should reflect the underlying pure state 
structure [HllEl. 

We have shown in particular that the multi-stage relaxation of the plaquette models 
leads to well-defined FDRs within each stage, with values ranging from X = 1 for quasi- 
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equilibrium regimes, non-trivial < X < 1 for local observables within activated stages, 
and non-trivial X < for global observables coupled to thermal activation. Therefore, 
plaquette models present an extremely rich aging behaviour that should be considered 
as an alternative paradigm to mean-field models when interpreting numerical data on 
physical systems such as disordered magnets or supercooled liquids. 
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Appendix A. Measuring response functions 

In this Appendix we give the generalisation of the 'no-field' method of |38[ IHHj to 
continuous time processes: we used this method to calculate all responses to local fields 
in this article. 

Suppose that we have a spin system with Hamiltonian H = Hq — ^ . hiSi, where 
Si G { — 1,1}; the sum is over all the spins, and the hi are time-dependent local fields. 
We choose a continuous time dynamics that obeys detailed balance, in which the spin 
on site i flips with a rate 

u/. = li ( A ]] 

' l + exp[/?(A, + 2/i,s,)]' ^ ■ ' 

where f3 is the inverse temperature; Aj is the change in Hq on flipping spin i; and 7, 
is a local rate that does not depend on the state of spin i or on the field hi. Glauber 
dynamics is the case when 7, = 1 throughout. In kinetically constrained models, we 
have 7j G {0, 1} according to state of the neighbours of spin i. In section ITTl we consider 
a modified Glauber dynamics, such that 7^ depends on the environment of each spin. 
However, the local field hi always couples to the dynamics as in ()A.1|) . 

The application the method of |HHl IHEj to this continuous time system is as follows: 
the integrated response of observable A to the field hi is 



{A{t)n,{t,t^))f,^=o (A.2) 



" d{f3hi 

where the field hi acts on site i between the times t„ and t, and 

where the sum is over the flips of spin i within the given time window; the notation t^^^ 
indicates that the summand is to be evaluated just before the spin flip. The strength of 
this method is that the small field hi does not enter (IA.3I). 
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We now outline a derivation of ()A.3|) . Let a trajectory be a sequence of 
configurations of the system at a series of time steps, where sequential configurations 
differ in at most one spin. The probability of a given trajectory is given by 



traj 



t/e 



(A.4) 



where represents the state of the whole system in the rth configuration; Po{Co) is the 
probability of the initial condition; e is the microscopic time step; and 



C'C 



1 



•i / i 



Wi 



(A.5) 



is the probability of a transition from state C to state C, in a time e. We use Si (or s[) 
to denote the state of spin i in configuration C (or C); the rates Wi are defined in ()A.1|) . 

We note that the microscopic time step must be small enough to ensure that Wcfi 
is always positive. Simple discrete time Monte Carlo dynamics have Wi < 1 for all z, 
and e = A^~^ where N is the total number of spins. Here we explicitly allow for e < A^~^ 
since we wish to take the continuous time limit. Also, recall that Ptraj is defined only 
for trajectories in which successive configurations differ in at most one spin. Thus we 
have Wc'c = 1, as long as the sum is taken over states C that differ from C in at 
most spin; this ensures conservation of probability. 

The definition of the stochastic average is then 



(^) — Araj -Ptraj, 

traj 



(A.6) 



where Atraj is the value of the observable A for the given trajectory, and the sum is over 
all possible trajectories. 

To calculate the integrated response, we simply write 



d{(3hi 



-Pt 



traj 



t/e 

E 

={twA)+i 

rl 

-W, 



t/e 

n 



Wt 



n n 

L-T ,L-T— 1 



d(/3/i, 



Po(Co), 



(A.7) 



where we assumed that the field acted only between times t^ and t, as before. We then 
write [diyc'c/d(/3/ii)] = Rcc^cc, which defines the matrix elements Rcc (at least for 
Wc'c > 0, which is the only relevant case). Hence, 

-(^) = y^Araj Ptraj ^ RCrfir-i- 

T=t„/e 



diPhi 



traj 



Explicitly constructing the matrix Rcrfi^-i as in Ref. 
continuous time, e — > 0, leads to (IA.3I). 



(A.8) 

and then taking the limit of 
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To calculate the defect response in the plaquette models, we also require the 
response for cases where the random field couples not to a single spin, but to a local 
function of several spins. For example, consider the response to a perturbation 

6H = -gi S1S2S3S4 = -giUi 

where the second equality defines the dual plaquette variable nj. In that case, 

-{Ait)) = {A{t)Si{t,t^)) (A.9) 



where the perturbation acts between times tw and t as before, and 

-2n. 



i.t„. 
■ flip 



1 + exp(-/?A„ ] 

aG{l, 2,3,4} ^ flipse(tw,t) "'*flip^ 



^^. 2n,.7a.exp(/jA. ^^^^^^ 



[l + exp(/3A„tO]2 

The sum inside the curly brackets is over flips of spin a, and the summand is again 
evaluated just before the spin flip. 



Appendix B. Field theory for defects in the SPM 

In this appendix, we define an effective reaction-diffusion system that takes account 
of the two growing length scales in the aging regime of the SPM. We calculate that 
correlation and response in this system at tree-level, following ^Hj- The methods are 
quite standard, so we quote only the main results; the review of |17] gives technical 
details of this formalism. 



Appendix B.l. Effective theory and dynamical action 

The effective theory is defined for defects on a cubic lattice of sites in d dimensions, 
with integer occupancy on each site. Diffusion takes place between neighbouring sites, 
and annihilation takes place on a single site. That is, 

UiUj -> (ni - 1), (nj + 1), rate ni Ai /g 
Ui — s> {ui — 2), rate nj(nj — l)Aj. 

The rates Dij and Aj depend on the local defect density. We also require their coupling 
to perturbations to the energy and to the inverse temperature (3. We write the energy 
as E = Ui — ^ ■ hiUi, and the appropriate rates are 



-I3(2~hi~hj) 



A. = 2Do (B-2) 

i' 

\ =Aoe-2^(i-'^»)5^^(>,. (B.3) 

i> 

Both rates are of activated form, with energy barriers close to two, since they require an 
intermediate stage containing a dimer. The factor in the denominator of 
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sets the coupling to the random potential h. The functions (7-,^] and g^^J decay with the 
distance between sites i and i', on a length scale i. They vanish aX i' = i and i' = j, 
and g^J^^j is symmetric in i and j. We also normalise such that gf}- = g^^^j = 1. 

The master equation for this system can be written in an operator formalism |48j . 
We introduce a bosonic algebra: [a^, aj] = [a], aj] = 0; [a^, aj] = and a vacuum state 
|0), such that aj|0) = 0, for all i. Let the probability of a configuration {nj} be P({nj}), 
and let |-P(t)) = J2{ni} ^{{''^i}^ ni('^l)"'|0)- Then the master equation is 

dt\P{t)) = -W\P{t)) (B.4) 

with 

W = 5^(l/2)(al - a])[(A, + D,,)ia. - a,) + (A, - A'Ol^^ + 

{ij) 

+ Y.^U{4f-M (B.5) 

i 

where Dij and Aj are now operators, since they depend on the densities uii = a\iaii. 

We now construct a path integral representation of the dynamics. We work in 
the coherent state representation, in which operators have matrix elements 0(0,0) = 
(0|e^«(^"+^)"'Oe^>'^'"'^|0)e~^>(^'+^)'^'. We take the continuum limit, denoting the lattice 
spacing by /q, and the position of site i by rj; we define 0(r) = Yli't'i^i.''^ " '^■d ^^id 
0(r) = {IqY 4'i^{.f~'fi)i so that the dimensions of the field 0(r, t) are those of density, 
and 0(r, t) is dimensionless. In a similar way, the field /i(r, t) = {IqY hi6{r — Vi) has 
dimensions of energy. The local density of excitations is 

p(r,t) = [l+0(r,t)]0(r,t). (B.6) 

We also define g{r) = gi'-i6[r — (r-j/ — r^)] and we choose gi'.ij to have the same spatial 
dependence: gi^-ij = g[ri' - (1/2) (r^ + rj)]. 

Finally, the path integral representation of the generating function for the dynamics 
is Z = 1 = / Vlcf), (j)] exp(— S'[0, (j)]), with 

S[^, 0] = jdt d'^r 0(r, t)dt^{r, t) + J dtH[cj), 0, t] (B.7) 

if [0, 0, t] = j d'^r d'^r' g{r - r') [1 + 0(r', t)]0(r', t) x 

|DV0(r, t) ■ [(1 + 2(3h{r, t)) V0(r, t) - 0(r, t)/5V/i(r, t)] 
+ A0(r, t) [2 + 0(r, t)] [1 + 2/3/i(r, t)]0(r, t)=^ 
+ 0(V^/^2)} (B.8) 

and A = Ao/g^ = A^o^'^ 

We work to linear order in the perturbation h and to quadratic order in a gradient 
expansion; we also neglect boundary contributions to the action. 

This concludes our definition of the effective theory. We now calculate the time 
dependence of the density, and the propagator in the system. 



where 
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Appendix B.2. Mean density and propagator 

For the calculations of this subsection, we set the perturbing field h = 0. This field only 
enters the calculation of the response, in the next subsection. To calculate the average 
density, we write 

dtn{t) = dt{p{r, t)) = {dt(P{r, t)). (B.9) 

where the average is over trajectories with initial density uq, and weights given by the 
action of (iRTll . 

Transforming to momentum space, we define (f)k{t) = J d'^r e*'°'^0(r, t), and 
similarly (f)k{t). The Fourier transform of g{r) depends only on k = \k\ and is denoted 
by Qk- Working at tree level, we evaluate {dt4>{r,t)) at the saddle point, whose position 
is given by the Euler-Lagrange equation 6S/6(j)^q{t) = 0. Writing out this equation for 
g = 0, we arrive at 

(9i0,=o(t)) = -A |d^rdV(7(r'-r)(A[0(t)]</)(r,t)20(r',t)) (B.IO) 

where y4[0(t)] is a function that depends on all of the Fourier components of (p, evaluated 
at a single time t. To evaluate the right hand side of ()B.in|l . we note that all expectation 
values of the form vanish (where again denotes a function that is 

local in time but depends on all Fourier components of the field </>). Further, the tree 
level approximation is that {B[(j){t)]) = for any function B. Hence, we arrive 

at 

dtn{t) = -2An(^)^ (B.ll) 

where we used the fact that gk=o is normalised to unity, and A[(j) = 0] = 2. Starting 
from an initial condition with density no, we arrive at the mean density: 

n{t) = ^° (B.12) 

We now calculate the propagator for the dynamics, defined by 

Gk{t,w) = {Nl^,)-\Mt)^~k{w)), (B.13) 

which satisfies lim^^^- w) = 1. To calculate the time dependence of Gk(t,w), 
we obtain the time derivative of (p-kit) from a second Euler-Lagrange equation, 
SS/Scpkit) = 0. At tree level, the only non- vanishing contributions to {(f)kit)dw4>-k{w)) 
are of the form 

{Mt)^-k{w)Uk[<P{w)]). (B.14) 

(Terms of the form (0fc(t)0_fc(w)0_fc/(w)[/^[0(w)]) vanish since every field must be 
contracted with a field at a later time.) Making the only possible contraction, and 
evaluating the function Uk[4>{w)] at the saddle point, we arrive at 

d^Gk{t,w) = Gk{t,w)Ukmw))] (B.15) 
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Finally, we use the Euler-Lagrange equation to identify Uk[{(p{w))] = Dk'^n{w) + 
2A(2 + gk)n{wY. Hence: 

w) = f e^'=(-)-^^W, (B.16) 

\n{w)J 

with 

Zkit) = Dk"^ /[2\n{t)\ (B.17) 

Armed with expressions for n{t) and Gkit, w) we are now in a position to calculate 
the correlation and response. 

Appendix B.3. Correlation and response 

The two-point defect correlation function analogous to Cd{k,t,t„) is 

C^{k, t, t„) = j d'^r e^'=-^[(p(r, t)p(0, t^)) - n{t)n{t^)]. (B.18) 

The integrand decomposes as 

e*'=-"[(0(r,t)0(O,tw)0(O,tw)) + (50(r,t)50(O,t„))], (B.19) 

where 6(j){r,t) = (f){r,t) — n{t). This correlation is to be evaluated a.t h = 0. 

At tree level, the integral of the first term of ()B.19|) is simply [Gk{t,t^)n{t„)]. 
We can evaluate the second term by direct consideration of its equations of motion. 
However, it is simpler to identify the tree level contributions to this quantity directly: 
expanding the time-ordered exponential e""^, the only terms contributing at tree- level 
come from: 

dsH[(f),^,s]y (B.20) 

This expectation value can be calculated straightforwardly, by making the allowed 
contractions and using the tree level relation {B[(f){s)]) = B[{(f){s))]. The result is that 



{6<p.kit)6<PkH) = -2A(Ar/^) / ds Gk{t, s)Gk{w, s)n{sf{l + 2gk) (B.21) 

Making the change of variables, Zg = Dk'^ /[2Xn{s)], leads to 

= - {Nl^oMw)Gk{t,w)zZ'-^^'' X 

/ dZs (1 + 2(7fc)22+23,g2.,-2.™ (g_22) 

This tree level calculation is valid in the aging regime, when n{t^) <^ n{t = 0). In this 
limit, the integrand is dominated by large Zg, and the behaviour is independent of the 
initial condition. We therefore evaluate the correlation at Zq = 0. 

The regimes of interest are /c£ 3> 1, for which ~ 0, and ki <^ 1 for which gk — 1. 
The integral can be performed in these limits. Bringing everything together, we arrive 
at 

G^{k, t, t^) = n{t„)Gk{t, t^)F[zk{t^),gk] (B.23) 
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where Zk{t) was defined in ()B.17|) . and 

F[z,0] = 1 - {1/A)[2z-^ - 2z-^ + ~ e-^')] (B.24) 

F[z, 1] = 1- (3/4) [2z~^ - Az'^ + Qz-^ - Gz'^ + 3^-^(1 - e'^')] (B.25) 

(Tliese functions are regular at z = 0, wliicli can be verified by direct expansion of tlie 
exponential.) 

The calculation of the response is also quite simple. We define the impulse response 

by 

6S 



h=0 

where /?.fc(tw) is the Fourier transform of the instantaneous perturbation h{r) which acts 
at time t^- 

Substituting directly for the derivative of 5*, the terms that contribute at tree level 
are of the form 

R^{k,t,t„) = {Nl^,)-\Mt)^-MWm.)]), (B.26) 

where (f)-k{t„)W[(f){t„)] contains the terms in ^^^^ ^ that are linear in 0. Evaluating 
this correlation function at the saddle point, we have i?</,(fc, t, tw) = Gk{t,t^)W[{(f){t^))]. 
and hence 

R^{k, t, t„) = 2X[zk{t^) - 2]n{t^fGk{t, t^) (B.27) 
Evaluating the FDR by taking a derivative of the correlation, we arrive at 

R<l>{ki t, ^w) 



^"^^'^'^^^ - d,^C,ik,t,t^) 

-2 + z,{t„) ^g_28) 



{l + gk + Zkit^)[l + -^]'jF[z,it^),gk] 



where we note that the FDR is independent of t, and depends on k and t^ only through 
the scaling variable Zk{t^). 

As discussed in the main text, the relevant limits for the SPM aging regime are 
H < 1, > p(tw), and < A;^ < p{t). These correspond to [zk{t) < l.gt = 0], 
[zk{t^) ^ 1,5'fc = 1], and [zk{t) <til,gk = 0] respectively. Calculating the FDR in these 
limits leads to 



R(f>{k, t, t^) 



X^(fc,t,M^ - = <^ l + 0{zk{t^)) [A;2»p(t„)] (B.29) 

' ' I _3_^0(^^(^^)_i^ [r^ « /t^ « p(t)] 

As expected for a reaction-diffusion system with activated dynamics [2Zj , the large 
k response has an FDR of unity, while the response at small k is negative. The 
intermediate regime with a well-defined negative FDR is a result of the presence of 
two length scales in the SPM, both of which grow with time. 
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